Interaction-induced stabilization of circular Rydberg atoms 
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We discuss a candidate solution for the controlled trapping and manipulation of two individual 
Rydberg atoms by means of a magnetic loffe-Pritchard trap that is superimposed by a constant 
electric field. In such a trap Rydberg atoms experience a permanent electric dipole moment that 
can be of the order of several hundred Debye. The interplay of electric dipolar repulsion and three 
dimensional magnetic confinement leads to a well controllable equilibrium configuration with tunable 
trap frequency and atomic distance. We thoroughly investigate the trapping potentials and analyze 
the interaction-induced stabilization of two such trapped Rydberg atoms. Possible limitations and 
collapse scenarios are discussed. 
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I. INTRODUCTION 

Among the many fascinating systems encountered in 
ultracold atomic and molecular physics are Rydberg 
atoms, i.e., highly excited atoms with large principal 
quantum number n. Their size can easily exceed that 
of ground state atoms by several orders of magnitude 
and at the same time is the origin of many extraordinary 
properties such as their massively enhanced response to 
external fields and, therewith, for their enormous polariz- 
ability |T]. In ultracold gases, the resulting strong dipole- 
dipole interaction among Rydberg atoms has been found 
to give rise to a non-linear excitation behavior: Ryd- 
berg atoms strongly inhibit excitation of their neighbors 
entailing a state dependent local excitation blockade l2)- 
[B], which on its part results in a collective excitation of 
many atoms [7H5]- Two recent experiments even demon- 
strated the blockade between two single atoms a few mi- 
crometers apart [11] . From an application oriented 
point of view, the strong dipole blockade effect renders 
Rydberg atoms promising candidates for quantum in- 
formation processing [T^ and allows the determination 
of the interaction potential of Rydberg atoms in a one- 
dimensional lattice . The large size of Rydberg atoms 
can also give rise to bonding interactions between Ryd- 
berg and ground state atoms. The scattering-induced 
attractive interaction binds the ground state atom to the 
Rydberg atom at a well-localized position within the Ry- 
dberg electron wave function and thereby yields giant ul- 
tra long-range molecules that can have internuclear sepa- 
rations of several thousand Bohr radii TT . The spectro- 
scopic characterization of such exotic molecular states. 
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named trilobite and butterfiy states on account of their 
particular electronic density, has succeeded recently [T3] 
and has triggered a revived theoretical and experimental 
activities [16,-.19j. 

Most of the experiments with Rydberg atoms still in- 
volve a large ensemble of atoms. They can therefore un- 
avoidably solely investigate effective and averaged prop- 
erties since individual atoms are typically not resolved. 
It is hence of great interest to study only a small number 
of Rydberg atoms that are preferably individually con- 
trollable and arrangeable with respect to one another. It 
is furthermore necessary to stabilize these Rydberg atom 
configurations against autoionization. An essential step 
in this direction is the trapping of electronically highly 
excited atoms. Several works have focused on trapping 
Rydberg atoms, based on electric [2D], optical [^TH^ . 
or magnetic fields [^^27) . Due to the high level den- 
sity and the strong spectral fluctuations with spatially 
varying fields, trapping or manipulation in general is a 
delicate task. This is particularly the case when both the 
center of mass and the internal motion are of quantum 
nature and the inhomogeneous external fields lead to an 
inherent coupling of these motions. 

In the present work we provide a candidate solution for 
the controlled trapping and manipulation of two individ- 
ual Rydberg atoms by means of a magneto-electric trap. 
Specifically, we consider a magnetic loffe-Pritchard trap 
that is superimposed by a constant electric field which 
induces a permanent electric dipole moment for the Ry- 
dberg atoms that can be of the order of several hun- 
dred Debye. As has been shown in a previous work, the 
resulting dipole-dipole interaction in conjunction with 
the tight radial confinement of the loffe-Pritchard trap 
gives rise to an effectively one-dimensional ultracold Ry- 
dberg gas with a macroscopic interparticle distance [28]. 
Here, we consider in addition the longitudinal confine- 
ment that arises for a non-Hclmholtz configuration of the 
loffe-Pritchard trap. In contrast to our previous work fo- 
cusing on the trapping of individual Rydberg atoms in 
two dimensions [251 - f?7] . this allows the controlled con- 
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finement of two single Rydberg atoms in three dimen- 
sions with variable trapping parameters and distance. 
We thoroughly investigate the resulting trapping poten- 
tials and analyze the interaction-induced stabilization of 
two such trapped Rydberg atoms. Possible limitations 
and collapse scenarios are discussed. 

In detail, we proceed as follows. In Section |ll] the 
Hamiltonian of a single Rydberg atom in the magnetic 
loffe-Pritchard trap superimposed by a constant electric 
field is derived and the corresponding adiabatic potential 
surfaces for the center of mass motion of the Rydberg 
atom are provided, as well as analytic expressions for the 
electric dipole moment of the Rydberg atom induced by 
the external electric field. In Section Hill we consider the 
interaction of two Rydberg atoms in the same trapping 
environment. In the regime of strong transversal con- 
finements, analytic expression for the equilibrium config- 
uration of the two atoms are derived. Weakening this 
restriction leads to three-dimensional potential surfaces 
and possible loss mechanisms that are investigated in 
the remainder of the section. Section Hvl outlines several 
routes to experimentally realize the proposed system. A 
brief summary is provided in Sec.|Vj In the Appendix we 
present a detailed derivation of the perturbative results 
for the induced dipole moment of the Rydberg atoms in 
the considered trap. 



external field, our initial Hamiltonian in the laboratory 
frame reads (atomic units are used except when stated 
otherwise) 
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The magnetic moments of the particles are connected to 
the electronic spin S and the nuclear spin S according 
to fjii = — S and — ~^vf~^' ^^^^ 9n being the nu- 
clear g-factor; because of the farge nuclear mass, the term 
involving ^2 is neglected in the following. 

The vector potential and the magnetic field of the lofi^e- 
Pritchard configuration read 
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II. THREE-DIMENSIONAL 
lOFFE-PRITCHARD CONFINEMENT FOR A 
RYDBERG ATOM 

A. Two-body Hamiltonian for a single alkali 
Rydberg atom 

In a highly anisotropic magnetic field configuration like 
that of a loffe-Pritchard trap, the strength of the mag- 
netic field can vary significantly over the extension of a 
Rydberg atom. The large size of Rydberg atoms can 
therefore modify the way they couple to the magnetic 
field compared to the coupling of ground state atoms. 
We incorporate the large extension of the atoms into our 
description by modeling a Rydberg atom by two parti- 
cles, namely, a valence electron (particle 1) and an ionic 
core (particle 2). This is particularly appropriate for al- 
kali atoms that are commonly used in Rydberg experi- 
ments. We include into our model the coupling of the 
electronic and the center of mass motion of the atom and 
hence do not resort to the infinitely heavy mass approx- 
imation. While the inclusion of the fine-structure and 
quantum defects can be readily done, it turns out not to 
be necessary for high angular momentum electronic state 
which we will be focusing on [29J. The coupling of the 
charged particles to the external magnetic field is intro- 
duced via the substitution, Pi ^ Pi — qiA{ri); qi is the 
charge of the i-th particle and A{x) is the vector poten- 
tial belonging to the magnetic field B{x). Including the 
coupling of the magnetic moments due to the spins to the 



The "traditional" macroscopic realization of the loffe- 
Pritchard trap uses four parallel current carrying loffe 
bars which generate the two-dimensional quadrupole field 
Bi that depends on the field gradient G. Encompassing 
Helmholtz coils create the additional constant field Be 
where B denotes the loffe field strength 30J. Bq des- 
ignates the quadratic term generated by the Helmholtz 
coils whose magnitude, compared to the first Helmholtz 
term, can be varied by changing the geometry of the trap. 
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R is the radius of the Helmholtz coils, and 2D is their dis- 
tance from each other. The geometry factor Q{D, R) van- 
ishes for 2D = R, which is known as the Helmholtz con- 
figuration. Exposing Rydberg atoms to a loffe-Pritchard 
trap in a Helmholtz configuration has been extensively 
studied in Ref. [21]. Here, we assume Q to be non- 
zero and positive, 2D > R. In this case, the absolute 
value of the magnetic field on the Z-axis, |B(0, 0,^)| = 
B\l + 2QZ^\, increases quadratically with \Z\. The ge- 
ometry factor reaches its maximal value, Qmax = To^^^' 
when 2D = y/^R, i.e., the smaller i?, the larger Q. 

Along the lines of Ref. [H] we employ the unitary 
transformation U — exp {| Be x r ■ il}, introduce rela- 
tive and center of mass coordinates (r = ri — r2 and R = 
(Miri + M2r2)/M with the total mass M ^ Mi+ M2), 
and omit the diamagnetic contributions. The Hamilto- 
nian describing the Rydberg atom in the loffe-Pritchard 
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trap becomes 

H^^^+HA + \Lr-B, + S- B{R + r) 

+Ai{R + r)-p + Aq{R + r)-p, (4) 

where we named the electronic hydrogenic Hamiltonian 
in field-free space Ha =p^/2 — 1/r. In Hamiltonian 
Q we neglect any contributions arising from the trans- 
formation of Aq{R + r) which is justified as long as 
\X\, \ Y\ <C As usual, the principal quantum number 
of the Rydberg state is denoted as n. 

For all relevant laboratory field strengths the spectrum 
of the Hamiltonian 0) is dominated by the field-free en- 
ergies Ea = — l/2n^hat are n^-fold degenerate. It has 
been shown in [55] that the inter-n-manifold couplings 
originating from the constant and the linear term in the 
magnetic field, Be and Bi, are negligible. The quadratic 
contribution fi\Bq{R) \ does not further constrain the pa- 
rameter regime where this approximation is valid. We 
can therefore restrict our study to a single sub-manifold 
with a given principal quantum number n. Considering 
the expressions 2{a'\xiPj\a) — eijk{a'\Lk\a) , we find 

H~^+t^-B{R)+n^+nQ + H' (5) 

where we named Hq :~ 2QZ{zLz — xSx — ySy + 2zSz), 
fx := \Lr -f and T-L-^ = G{xypz + xSx — ySy) as in 
Ref. '29]. The terms 

H' = S- Bq{r) - 2Q{zXSx + zYSy + {xX + yY)Sz) 

+ j[{x'^y + y^- 4yz^)px + {-xy^ - x^ + 4xz^)py 

+ {-y'^Py - "ix^Py + ^z^Py + 2xypx)X 

+ ( x'^px + iy'^Px - 'iz^Px - 2xypy)Y] (6) 

are small corrections as long as 4^n^ <C 1. For 
B = lOG, G = lTm~^ and n = 30 this condition reads 
Q <^ 1.5 X 10~^^. To reach geometric parameters Q as 
large as 10~^^, the coils of the loffe-Pritchard trap would 
have to be as close as 12 fj,m. For a macroscopic trap, 
the above condition is therefore always valid. In con- 
trast to the finite-size term "Hj, the term Jig depends on 
the center of mass position. Since Lz and Si are diag- 
onal in the hydrogen basis, the latter is proportional to 
a dipole matrix element and hence has only off-diagonal 
matrix elements. Comparing its second order energetic 
contribution with IQL^Z^j (which is a part of fi ■ B), 
and assuming the energetic gap of adjacent surfaces to 
be B/2, cf. Ref. [2^, we find it to be negligible as soon 
as An^Q ^1. In macroscopic traps this restriction can 
only be broken with principal quantum numbers n of the 
order of several hundreds. 

Our working Hamiltonian thus reads 

niP=^ + ti-BiR) + n^=:^ + n.. (7) 



In order to solve the remaining coupled Schrodinger equa- 
tion, we adiabatically separate the relative and the center 
of mass dynamics by projecting Eq. ([7| on the electronic 
eigenfunctions ip^ that parametrically depend on the cen- 
ter of mass coordinates: 

n,\v^{r;R))=E^{R)\^,{r;R)). (8) 

We are thereby led to a set of decoupled differential 
equations governing the adiabatic center of mass motion 
within the individual three-dimensional energy surfaces 
E^{R), i.e., the surfaces Ei^{R) serve as potentials for the 
center of mass motion of the atom. The non-adiabatic 
(off-diagonal) coupling terms that arise within this pro- 
cedure in the kinetic energy term can be neglected in 
our parameter regime since they are suppressed by the 
splitting between adjacent energy surfaces |29j . 

B. Electronic potential energy surfaces 

Approximate expressions for the potential energy sur- 
faces Ef^{R) can be found analytically when the ratio of 
the magnetic field gradient and the loffe field is small, 
more specifically if n?G/B <C 1. In this case, the finite 
size term T-L-^ is negligible compared to the contribution 
of /i, • B{R). The latter can be diagonalized by applying 
the spatially dependent unitary transformation 

where tana = By{Bl + B^)-^^^ and tan^ = -B^/Bz- 
The spatial dependence is introduced by evaluating the 
magnetic field components Bi at the Rydberg atom's cen- 
ter of mass position, i.e., Bi = Bi{R). We note that the 
field-free Hamiltonian Ha is invariant under the trans- 
formation U . Moreover, 

U^lU^B=]^{Lz + 2Sz)\B\, (10) 

where Lz and Sz are now defined with respect to the 
local quantization axis. This procedure is equivalent to 
rotating the system into the local magnetic field direc- 
tion. The adiabatic potential energy surfaces hence read 

E,{R)^C^+ms)\B{R)\. (11) 

Expanding the absolute value of the magnetic field 
around its minimum in the trap center, 

|B(0, 0, Z)| w 2QZ^ ^ B{1 + 2QZ^) , 

\B{X, YM^B+{^-Q^p' + 0{p') , (12) 

yields the harmonic confinement known from ground 
state atoms in a loffe-Pritchard trap [p = y/X"^ + F^) . 

In the considered limit, n?G/B — )■ 0, the energetically 
uppermost electronic adiabatic potential energy surface 
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is the only non-degenerate one. The electronic state that 
corresponds to the uppermost surface is (in the rotated 
frame of reference) the circular one whose angular mo- 
mentum projection quantum number mi — I — n — 1 is 
maximal within the given n-manifold. Hence, the trap 
frequency experienced in this surface exceeds the one of 
a ground state atom by a factor n — 1. This entails ex- 
tremely large transversal trap frequencies for the external 
motion such that the extension of the center of mass wave 
function can become even smaller than the extension of 
the electronic cloud of the Rydberg atom [53]. Since the 
uppermost surface additionally suffers from the smallest 
non-adiabatic couplings and due to its non-degeneracy, 
it is best suited for a controlled confinement. 



In the derivation of Eq. (11) we neglected the finite 



size term Ti^ since it only involves relative coordinates 
and therefore constitutes to first order solely a constant 
energy offset to the surfaces. Because of the coupling of 
the relative and the center of mass motion, however, its 
contribution to the electronic energy and the electronic 
wave function will ultimately depend on the center of 
mass coordinates. The resultant admixture of other hy- 
drogenic states to the circular state entails a non-zero 
permanent electric dipole moment outside the trap cen- 
ter, which will be discussed below. The corresponding 
energy surface itself, on the other hand, hardly shows 
any deformation even for large magnetic field gradients. 



C. Electric dipole moments 



Let us proceed by studying the expectation value of 
the electric dipole moment of the Rydberg state in the 
uppermost adiabatic energy surface. A non-zero dipole 
moment arises from parity symmetry breaking terms in 
the Hamiltonian. The only such term in the working 
Hamiltonian (It]) is the finite-size term Ti^, whose im- 
plicit i2-dcpcnacnce will eventually entail a spatially de- 
pendent dipole moment. In the following, we pursue a 
perturbative treatment of Hj that gives rise to an ex- 
plicit expression of the resulting dipole moment. To this 
end, the off-diagonal matrix elements of the perturbation 
operator ~ G need to be much smaller than the corre- 
sponding unperturbed energy level spacings AE ^ \B\. 
This yields the requirement B/G 3> n'^ which is easily 
satisfied for typical loffe field strengths. 

The perturbative treatment is detailed in Appendix [X] 
It results in the expression for the permanent electric 
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where x := 9n2 (2n2 - 3n - ^4^2 _ iQn 6 + 1)/(8A£;) , 
A — G/3, and AE « \B\/2 is the energetic gap between 
the uppermost surfaces at the trap center. As can be de- 
duced from Eq. (13), the electric dipole moment is per- 



pendicular to the local direction of the magnetic field, 
■ B — {). It vanishes on the Z-axis. 
Further control of the electric dipole moment, both re- 
garding the magnitude as well as the steric properties, 
can be gained by applying an additional electric field. In 
the following, we thus consider a modified loffe-Pritchard 
trap with an additional electric field F — (-F!r, 0, 0) point- 
ing in the x-direction as in Ref. [5S]. The latter can be 
treated perturbatively as long as ^ B/n. The elab- 
oration of the perturbative treatment (presented in Ap- 
pendix |A| shows that the energetic contribution of the 
electric field Hamiltonian is of second order, 

A2,e(2.F.) ^ ^n\n - l){cos fi'' + sin fi''), (14) 

where Xp — |F|/|B|. For vanishing Q and with the 
approximate expression for the energetic separation be- 
tween the coupling surfaces, AE sa \B\/2, this simplifies 
to 

-^Kn[n 1)^2 ^ ^2^2 + G2X2 • ^^^^ 

The perturbative contribution to the uppermost surface 
due to an external electric field is thus positive and it is 
maximal on the Z axis. For small atomic displacements 
\R\ ^ 5/G, it can be considered a mere offset to the 
uppermost surface. 

The electric field induces an electric dipole moment 
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F = -—n (n-l)^ 
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'B2 



5^ + 52 

— B^By 
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that depends on the ratio of the field strengths, as ex- 
pected, and on the cubed principal quantum number. 
Surprisingly, however, only on the Z-axis it points along 
the direction of the generating electric field. Similar to 
d^, dp is in general perpendicular to the local quanti- 
zation axis, which is set by the magnetic field direction, 
i.e., dp ■ B ~ d^ ■ B — 0. In addition, d^ and dp are 
perpendicular to each other on the Y-axis; on the posi- 
tive A"-axis they are parallel while being anti-parallel on 
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the negative X-axis. These properties only apply as long 
as the perturbative treatment is applicable. Our numer- 
ical results show that the dipole moment aligns with the 
electric field for a larger electric field strength. 

Analyzing the symmetry properties of the combined 
electric dipole moment d = + dp reveals that even in 
case of a longitudinal confinement an approximate sym- 
metry in Z remains, d{X,Y,Z) sa -d{X,Y,-Z). Much 
stronger than the .Z-dependence of the electric dipole mo- 
ment is its dependence on the transversal coordinates X 
and Y. If no electric field is present, the relevant electric 
dipole moment is d^, generated by the finite size term 
H-y. The symmetry properties of its components in the 
XY-pl&ne read 



-^72 



{X,Y) 




-X,-Y) 



(17) 



We have furthermore d~f,z{X,Y) = — d-y ^(±Y, ±X) and 
dj^xiX,Y) = d^^y{Y, X). The dipole moment induced by 
the external electric field, dp, exhibits different symme- 
tries, 





-X,-Y) 



(18) 



The sum of the contributions is therefore only symmetric 
with respect to a reflection about the X-axis {Y — ?> —Y)- 



d=\dy\ {X,Y) = 




iX,-Y) 



(19) 



This can be seen in Fig. [T] where the components and 
the absolute value of d are depicted. The parameters 
are chosen such that d^ and dp are of the same order 
of magnitude. Already for moderate electric fields, how- 
ever, dry is a mere perturbation to d and the symme- 
try properties of d are approximately those of dp, cf. 
Eq. ( 18 ). We note that the correct procedure in the com- 



bined magneto-electric trap is to consider the perturba- 
tion operator consisting of the sum of the electric field 
term Hp and the finite size term "H-y, rather than adding 
up the dipole moments d^ and dp generated by the in- 
dividual contributions. However, as elucidated in Ap- 
pendix |Xj the latter approach is exact in first order. We 
found very good agreement of the perturbatively com- 
puted expectation values for d with the numerically cal- 
culated values, both for vanishing as well as for (small) 
finite electric field strengths. 



III. INTERACTION-INDUCED 
STABILIZATION 

In this section, we extend our studies by considering 
two Rydberg atoms A and B that are trapped in a loffe- 
Pritchard trap and that interact via their electric dipole 
moments. The adiabatic Schrodinger equation for the 
two-atom center of mass wave function, j^'y^s), for this 
situation reads 

[Ta + Tb + Va{Ra) + Vb{Rb) 

+ ^„t(i2A, Re)] \^ab) = E\^>ab) , (20) 

where VSnt contains the interaction energy that depends 
on the positions of both atoms. The one-atom potential 
for the atoms A and B, Vyi(R) — Vb(J^), can be approx- 
imated for high-Ioffe-configurations by the analytically 
diagonalized term (10). The interaction potential l^nt 
will be discussed in detail in subsection IIII Al Because 
of the strong transversal confinement in the considered 
loffe-Pritchard trap, we restrict our considerations in a 
first step to the Z-axis in subsection |III B| In this simpli- 
fied geometry we analytically find a stable configuration 
of the atoms in which their distance is easily tunable 
without affecting neither stability nor trap frequencies. 
In subsection IIII Cl we extend our considerations to three 
dimensions and dwell on the question of stability. The 
last subsection is dedicated to experimental implementa- 
tions suggesting different ways of realizing the system. 



A. Rydberg-Rydberg interaction 

The interaction energy T4nt of two Rydberg atoms - 
each modeled by a core and an electron - can be formu- 
lated using the electric dipole moments of the individual 
atoms as long as the inter-atomic distance is large com- 
pared to the distance of the electrons to their respective 
cores. To this end, we write the Coulomb interaction 
between the charges of the different atoms. 



V\nt{rA,rB,RAB) _ 1_ 

e2/4^eo " |fi 



1 



AB\ 



\Rab - fsl 

1 



I-Ras+^aI \RAB-{rB-rA) 



(21) 



as a multipole expansion in the small parameter Ajnt = 
{ta.b) /Rab- Here, we abbreviated the vector connecting 
the ionic cores by Rab Ra—Rb and Va/b denotes the 
electronic relative vectors with respect to the cores A, B. 
Due to the neutrality of the interacting constituents, the 
only non-vanishing term up to third order in the expan- 
sion of T^nt is the dipole-dipole term Vdd- If we abbre- 
viate the projections of the electronic coordinates onto 



the vector connecting the cores as 



' AB ■- 



{rA -rB) ■ Rab 



''B' ^ab 



i ■ Rab and 

Rab /Rab, 



the multipole terms up to fourth order in the expansion 
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The only symmetry that survives when a small electric field is present is the symmetry with respect to the X-axis, 
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Parameters: B — 0.1 G, G 
atomic units eao = 2.54 Debye. 
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a.u.= 5.14 • 10 Vcm . Electric dipole moments are given in 



of the interaction potential, 

Vintir A, rB,RAB) = Vm + Klq + 0(A^ 



can be rewritten as follows: 
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(23) 



The last line in ( 23 ) holds if rj^~r^, e.g., for two circular 
Rydberg atoms in the same n-manifold. In this case the 
dipole-quadrupole interaction Vdq vanishes if van- 
ishes, that is when the electric dipole moment expecta- 
tion values for both atoms are identical. 

If the interaction operator T^nt is treated as a pertur- 
bation to the electronic Hamiltonians of the individual 
Rydberg atoms, Ha and Hb, it is favorable to repre- 
sent Vint in single-atom electronic eigenstates. We hence 
use the two-electron basis {\ff',ff)} = where i 

and j number the single-atom adiabatic electronic wave 
functions in the rotated frame of reference: 

{Ha + HBM;^f) - (E, + E,)\^t;^f) . (24) 

Note that we omitted the antisymmetrization of the two 
electrons, which is valid in the asymptotic region we are 
considering where the electrons are well localized at the 
respective Rydberg atoms. The leading order of Vint is 
given by the dipole-dipole interaction operator which can 
be represented in the above basis as 



RlB{^'f\VM\tj) 



= {^'\r\^)■{f\r\J)-3{^'\r^\^){f\r^\J) 



di>i{RA) ■ dj'j{RB) — 3d^i{RA, RAB)dff j{RB: Rab) , 

(25) 



where df^ 



= (ilrlj) ■ Rab- For configurations 
close to the Z-axis, i.e., when Rab ~ (0,0,1), the 



last term can be approximately written involving the 
z-components of the electric dipole moments only, 

di^iiRA, RAB)dj',j{RB, Rab) ~ di'i^z{RA)dj'j^z{RB)- 

For finite interaction between the atoms, the two-atom 
basis states are no longer eigenstates of the sys- 
tem. Looking at the state where both atoms are circular, 
l^cSV'c) = the only non-vanishing transition dipole 
matrix elements are 



(ll|14d|33) 



9n2(n- 1) 



Rab 



(26) 



where |3) denotes the state with I — mi — n ~ 2 in the 
rotated frame of reference. This coupling is small as long 
as 



(II|ydd|33) <^5 ^ Rab'^ti 
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1/3 



= Rc (27) 



where 5 ~ \B\ is the energetic separation of the surfaces. 
For the parameters i? = 10 G and n = 30 this yields 



Rab 2> 1.3 /im, which allows us to use the form (25) in 
first order in the following, i.e., only considering diagonal 
elements. 



B. One-dimensional stable configuration 

A loffe-Pritchard trap can provide an extremely strong 
confinement for Rydberg atoms in the transversal, i.e., 
XF-direction [SS]. We now want to take advantage of 
this peculiarity in order to restrict the study of the total 
potential T4ot '■— Va + Vb + Vdd to the Z-axis. In Section 
|III C| we investigate the requirements on the magnetic 
field parameters to guarantee that this simplification is 
permitted. There we find that for large enough gradi- 
ents G this is always the case since they entail strong 
transversal confinement. To simplify the situation even 
further we impose an external electric field pointing in 
the X-direction that keeps the atoms away from each 
other and prevents autoionization [28]. As before, the 
following discussion focuses on the uppermost potential 
surface, emanating from the circular state. 
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Small oscillations of generalized coordinates 

In order to study the one-dimensional configuration we 
set the coordinates X and Y to zero and assume Q to be 
non-zero and positive which generates the confining po- 
tential in the Z-direction. Please note that in this case 
rf = Ti ■ Rab = Zi in Eq. (23), since the atomic sep- 
aration vector coincides with the Z-axis. We calculate 
the expectation value of the electric dipole moment of 
the single-atom eigenstate |V'c) via Eq. (25). Addition- 
ally inserting the trapping potential Eq. ( 12 ), naming the 
atoms such that > Zg, and omitting the constant po- 
tential offset 2nB, the two-atom potential represented in 
the single-atom electronic eigenfunctions reads (again in 
the rotated frame of reference; 47reo/e^ 



Vr 



'{Za,Zb 
--2nQ{Z\ 



1 a.u.) 



+ zl) 



' 4{B + 2QZI){\Za- ZsiriB + 2QZI)- ^^^^ 

Utilizing generalized coordinates for the distance of the 
atoms and for their center of mass, Z^ = Za — Zb > 
and Zs — {Za + Zb)/2, respectively, the total potential 
( 28 ) translates to 



iQ {Z 



l + ^zl) 



vt\Zo.Zs)^^-^n\n~lf^^^^ 

(29) 

Here we approximated B + 2QZ\ g ^ B, which is valid 
as long as 



\Zs\+Zd/2^ VB/{2Q) . 



(30) 



The first term in Eq. (29) is the approximate version 



of the dipole-dipole interaction operator. It only depends 
on the distance of the atoms. Higher order terms orig- 
inate in the quadratic Z-dependence of the interacting 
electric dipole moments. They become significant only 
for very large Z or exceptionally strong parameters Q 
reachable on atoms chips. The coordinate for the cen- 
ter of mass of both atoms, Z5, appears as the quadratic 
shift AnQZg. An equilibrium configuration of the atoms 
is therefore bound to be symmetric around the origin, i.e., 
Zs = 0. Minimizing the energy of the two-atom poten- 
tial within this approximation, we find the equilibrium 
position at 



, 



z 



D.min 



F2 



8B^Q 



23/5 y B^Q 



(31) 



The expression for the equilibrium distance, Z/^ min, can 
hence be readily controlled by the electric field strength 
Fx ■ The condition of validity of our approximation ( 30 ) 
at the equilibrium position (31 1 reads 



CD 



1 



2B^ 



3n 
T 



2B9 



< 1 . 

(32) 



TABLE I. Explicit values for cd, Eq. (321, which measures 



the qu ality of the approximation ( 30 1 at the equilibrium po- 
sition (31 1. The values are computed for the geometry pa- 



rameter Q = 6 X 10~ which is around the highest values 
reachable with macroscopic lofFe-Pritchard traps [STl . The re- 
striction CD <C 1 is violated only for very low loffe fields paired 



with 


electric 


fields that would 


ionize the Rydb 


erg atoms. 


(10- 


" a.u.= 


5.142 206 32 V/cm) 






F 




10-1^ a.u. 


10"" a.u. 


10-1" a.u. 


B = 


0.1 G 


0.0297 


0.0747 


0.1877 


B = 


IG 


0.0075 


0.0188 


0.0471 


B = 


10 G 


0.0019 


0.0047 


0.0118 


B = 


100 G 


0.0005 


0.0012 


0.0030 



See Tab. |T] for explicit values. 

Since we are interested in the motion of the system 
around a stable equilibrium configuration, we expand 
the potential in a Taylor series around that equilibrium 
and solve the corresponding classical eigenvalue problem. 
The resulting frequency for the center of mass and rela- 
tive motion read in the harmonic approximation 



20nQ 



4nQ 



(33) 



2M ' ^ 2M ' 
where 2M is the total mass of the system. It is worth 



noting that within the approximation ( 30 ) the eigenfre- 
quencies are independent of the electric field strength F^ . 
They indirectly depend on the loffe field strength since 
Q = BQ. 



Tuning the distance of the atoms 
The fact that the equilibrium distance Zu , 

Eq. m, 



of the 



atoms 



can be increased without changing 
the trap frequency by just increasing the electric field 
strength F is depicted in Fig. |2][a) . The two-atom poten- 
tial and its harmonic expansion around the equilibrium 
position are drawn for different electric field strengths. 
In subfigure [2jb) a magnified view of the minimum is 
provided, demonstrating the validity of the harmonic ap- 
proximation. The figure shows that we can safely assume 
the center of mass ground state to be a Gaussian with 
the corresponding trap frequency cobi- Subfigure [2][c) 
in addition shows expansion coefficients of the potential 
around the equilibrium position (computed without as- 
suming the approximation ( 32 ) to be valid) . It is evi- 



dent from the plot that the harmonic description of the 
potential around the equilibrium position is a good ap- 
proximation for small loffe field strengths. Therefore, 
not only the trap frequency but also the center of mass 
ground state remains the same for different values of Zb . 
To give a numerical example, the variation of the trap 
frequency lub does not exceed lO"*^ for a loffe field of 1 
Gauss as long as the electric field strength is smaller than 
F« 2.3 V/cm. 
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FIG. 2. (color online) (a) Two-atom potential energy, 
Eq. | |28[ ) (solid), the harmonic approximation (dotted) and 
the ID ground state energy along Zd for Zs = and dif- 
ferent electric field strengths F. Parameters: B = 10 G, 
g = 6 X lO'^'^B, n = 30, F = 0.0514 V/cm (black), 
0.257 V/cm (blue/gray), 0.514 V/cm (yellow/light gray), (b, 
c) Quality of the harmonic approximation for the first set 
of parameters. Subfigure (b) depicts the harmonic approxi- 
mation (dotted) of the two-atom potential Viot.circ (solid) in 
a region around the local minimum (filled: Gaussian). The 
double-logarithmic plot in subfigure (c) shows the quadratic 
(dotted), cubic (dashed) and the quartic (dot-dashed) coeffi- 
cient of the expansion of the potential around the equilibrium 
position as a function of the loffe field strength. 



Quadrupole-quadrupole repulsion 

In this subsection we study the influence of higher or- 
der niultipole interactions and answer the question in 
which situations they can change the behavior of the sys- 
tem. We do this on the Z-axis for vanishing electric field 
and consider - as throughout the present work ~ the cir- 
cular Rydberg state. 

When no electric field is present, circular Rydberg 
atoms located on the Z-axis of a loffe-Pritchard trap do 
not exhibit a permanent electric dipole moment. They 
are hence neither subject to dipole-dipole interaction nor 
to dipole-quadrupole interaction. Since the circular elec- 
tronic wave function is not spherically symmetric, they 
feature a quadrupole moment, however. This is a first- 
order effect and the repulsive quadrupole-quadrupole in- 
teraction is hence the leading order of the interaction 
potential. For the special configuration we are consider- 
ing, it can be calculated using the simplified expression 



qq 



4Z 



22 r('22,22N -|t;22 



D 



+ '2{xaXb + VaVb - 4:ZAZBf 



(34) 



With the matrix elements (V'ck^lV'c 



i^Pclx^l^Pc) = {^c\y^m = in3(n+ I), and (V'ck'IV'c) = 
{n + 1), we find 



1^2 



^qq=^^^*(^+l)'(^+^)' 



1 3 
^2' 



(35) 



For low enough center of mass kinetic energy, the re- 
pulsion of the atoms due to the quadrupole-quadrupole 
interaction could in principle stabilize Rydberg atoms on 
the Z-axis against auto- ionization. We must not forget, 
however, that the van-der-Waals coupling as a second 
order contribution to the multipole interaction can be of 
similar strength for low enough distances in the consid- 
ered parameter regime. 

The situation changes completely when an electric field 
is applied. The induced dipole moments scale linearly 
with the field strength and the dipole-dipole interaction 
then depends quadratically on F. From the first term in 



Eq. ( 29 ) its magnitude can be estimated to be 



dd,2-axis 



1 /9 3F 
zh V2 B 



(36) 



The dipole-quadrupole interaction happens to be zero on 
the Z-axis even for finite electric field strength. In or- 
der for the dipole-dipole interaction, |Vdd| ^ I'iP/^iii to 
dominate the quadrupole-quadrupole interaction, \Vqq\ ^ 



^n^/Z^, the condition 







fZDF\ 






y nB J 



(37) 



must be fulfilled. For a loffe field strength _B = 10 G and 
an electric field strength as low as F = 10^^^ a.u. (and 
n = 30) this reads Zjj ^ 35000 a.u.= 1.7/im. Increas- 
ing the electric field strength to F = 2 x 10^^^ a.u. al- 
ready yields Zd ^ 1700 a.u.= 90 nm. For the examples 
above it is therefore legitimate to neglect the quadrupole- 
quadrupole interaction. Around the equilibrium config- 
uration Zu ^in of the atoms, Eq. (31), condition (371 is 
even easier to fulfill. 



C. Three-dimensional stable configuration and 
collapse 

Very strong transversal confinement leads to the one- 
dimcnsional situation discussed in the preceding Chap- 



ter III B Since the magnetic field gradient G only influ- 
ences the transversal but not the longitudinal confine- 
ment, the respective trap frequencies can be altered in- 
dependently. If the transversal confinement is decreased 
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FIG. 3. (color online) Sections through the six-dimensional two-atom potential around the local minimum. The Zo-coordinate 
of the minimum is indicated by the vertical lines. The thick black contours are plotted at the energy of the two-atom harmonic 
center of mass ground state which corresponds to half the sum of all six trap frequencies, Eq — 26 kHz. The contour plots are 
chpped at energies 125 kHz higher than the minimal energy configuration. The dashed and dotted lines in the plots on the left 
hand side indicate the quality of the single-atom-surfaces-approximation introduced in Sec. |IH A| They are drawn where the 
ratio of dipole-dipole interaction energy and energetic distance of adjacent single-atom electronic surfaces, Edd/^E, equals 0.1 
(dotted) and 0.01 (dashed). Parameters: B = 30 G, G = 10Tm-\ Q = 6 x 10""B, F = 2 x 10~" a.u.= 10.28 V/m, n = 30. 



and/or if the longitudinal confinement is increased the 
local minimum of the potential mentioned above turns 
into a saddle point: The tendency of the dipole-dipole 
interaction to force the atoms to step out of the Z-axis 
wins against the confining nature of the transversal mag- 
netic field gradient. The atoms then attract each other 
and most probably eventually ionize. These statements 
are substantiated and refined in the following. 

The two-atom interaction potential exhibits an approx- 
imate longitudinal symmetry as can be justified on ac- 
count of the smallness of Q = Q/B. More precisely, 



the conditions 2QZ'^ < 1 [Eq. (30)] and 4Q|Z| < G/B 



have to be met. Then, the dipole-dipole interaction Vdd 
between two atoms in the circular state exhibits a depen- 
dency on Zs that is negligibly small. The confinement 
due to the lofFe-Pritchard magnetic field configuration, 
on the other hand, is harmonic around the origin in all 
directions. Within the range of validity of the above ap- 
proximations, we can thus conclude that minimizing the 
total energy of the two atoms always leads to a symmet- 
ric configuration, Zs = 0. We therefore set Zs — for 
the following analysis. 

In order to characterize the six-dimensional adiabatic 
two-atom potential around the local stable minimum we 
first keep a strong transversal confinement which is quan- 
tified by the trap frequencies in X- and K-direction. For 
a single atom in a loffe-Pritchard trap and for finite Q 
they read 



UJ 



X 



2n - 2Q{B + 2GZ) 
M B + 2QZ2 
2n G^ - 2Q{B - 2GZ) 
M B + 2QZ2 



(38) 



Corrections due to the dipole-dipole interaction energy 
when two atoms are considered are proportional to 



,-5/2 



for a dominating loffe field. 



M B^D 

In what follows, let us investigate the exemplary pa- 
rameter set B = 30G, G = 10 Tm'^, Q = 6 x IQ-^^B, 
F = 2x 10^^^ a.u., and n — 30. All the inequalities, that 
have been formulated up to now in order to measure the 
quality of the applied approximations, hold with a con- 
fidence factor of at least 10^. The only exception is the 
quality of the perturbative approach to determine the 
electric dipole moment. The corresponding requirement 
explicitly reads nF^/B = 0.047 <C 1 for the chosen pa- 
rameters which is still satisfactory. 

Diagonalizing the Hesse matrix [d'^Vtot /{d Rid Rj)] at 
the local minimum position of the total, i.e., two-atom 
potential and extracting the trap frequencies along the 
principal axes from the eigenvalues and its eigenvectors, 
respectively, yields 



ll.lkHz, 
ll.lkHz, 
10.7kHz, 
9.7kHz, 
5.0kHz, 
4.5kHz, 



{Xa ~ Xb ) , 
{Ya = -Yb), 
{Xa = -Xb) 
{Zd), 



(39) 



where the equations in brackets define the directions of 
the principal axis. It is convenient to introduce appro- 
priate generalized coordinates for all three spatial dimen- 
sions, namely, Rd — Ra — Rb and Rs — \{Ra + Rb)- 
Sections of the total potential around its minimum are 
shown in Fig. [3] 
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The trap frequencies in (391 have to be compared to 



the radiative lifetime of the interacting atoms. The field- 
free lifetime of the electronic state corresponding to the 

uppermost surface is r(n,n — 1) ~ jc^ ia)^ ~ 2.3ms 
[55] , Even the reduction of the lifetime due to admix- 
tures to the pure circular state originating from the fi- 
nite size term and due to the coupling to the electric 
field leaves it close to the field-free value j28J. We hence 
expect more than 10 oscillations of the atomic motion to 
be observable within the lifetime of the Rydberg state. 
Since ~ nQ = nBQ the trap frequencies a; can ob- 
viously be increased by the choice of a larger principal 
quantum number n, with higher loffe field strengths B 
and by allowing for a stronger longitudinal confinement. 
The latter can be achieved by shrinking the trap onto 
an atom chip [551. We emphasize that the stable two- 
atom configuration is not immediately lost when one of 
the atoms decays to the circular state of the adjacent 
n-manifold since the electronic properties of that state 
are very similar and, therefore, so are the electric dipole 
moment and the adiabatic surface. 




FIG. 4. (color online) Spatial dependence of the direction of 
the electric dipole moments and dp (for Q — 0). (a) The 
electric dipole moment that originates from the finite size 
of the Rydberg atom. It is perpendicular to the local magnetic 
magnetic field direction (yellow arrows) and it vanishes on the 
Z-axis. On the positive (negative) X-axis it is parallel (anti- 
parallel) to the electric dipole moment dp that is induced 
by the electric field and depicted in the subfigure (b). (c) 
Orientation of dp with respect to the vector Rab connecting 
the interacting atoms A and B when they are displaced in 
the same X-direction, Xs 7^ (solid line), or in different X- 
directions, Xd 7^ 0, (dashed lines). 



Loss of confinement and collapse 

We now study mechanisms that endanger the stabil- 
ity of the equilibrium position on the Z-axis when the 
transversal confinement is relaxed. We identify two sit- 
uations in which this happens. One of them is the loss 
of the confining property of the loffe-Pritchard field con- 
figuration for a single atom for large ratios B/G. Fur- 
thermore, the stability of the equilibrium configuration is 
also lost as soon as the transversal confinement becomes 
smaller than the transversal anti-confincmcnt due to the 
dipole-dipole interaction. 

Regarding the first mechanism, the transversal mag- 
netic confinement for a single atom is only guaranteed as 
long as the respective curvature is positive. This yields 
the restriction 



52 



> 20 (1 



-4-|Z| 



(40) 



In order for this condition to be broken at the origin of 
the trap, the ratio G/B must be extremely small since 
the highest reachable values for the geometry parameter 
Q in macroscopic loffe-Pritchard traps are around 10~^^ 
(we use Q = 6 X 10^^^ for all presented examples [H]). 
For large enough displacements in Z-direction, however, 
the condition can always be broken. To give a sense of 
the numbers, we insert the exemplary parameter set B = 
lOG, G = 2Tm~i and g = 6 x IQ-i*^ to find that the 
displacement \Z\ must be as large 2 x 10^ a.u. — 1 mm to 
break the condition (40). We do not consider such large 



atomic distances from the trap center for any example. 

The second reason for the loss of the stable equilibrium 
configuration on the Z-axis is the dipole-dipole interac- 
tion between the two atoms. Besides being the interac- 
tion of longest range between neutral atoms, the major 



property of the dipole-dipole interaction is its anisotropic 
character. This comes into play when the atoms can 
step out of the Z-axis and the angles between the elec- 
tric dipole moments and the connecting vector change 
(Fig.|4|). 

Polarized case: In order to simplify our considera- 
tions let us first assume that the external electric field 
F ~ {Fx, 0,0) fully polarizes the atoms. This happens 
for relatively large electric field strengths as discussed 
in App. |A 2[ Then all dipole moments point in the X- 
direction. For the sake of clarity we now additionally 
assume that the magnitude of the dipole moments does 
not depend on the position of the atom in the trap. Then 
a displacement of both atoms A and B from the Z-axis 
in the same direction, Xs ^ and/or Ys ^ 0, does not 
change the interaction energy since the angle 9 between 
the two dipoles does not change. The angle also stays 
the same for a displacement of the atoms in opposite Y- 
directions. Yd 7^ 0. Here the interaction energy decreases 
only slightly due to the increase of the distance Rab 
of the atoms. A displacement of the atoms in opposite 
AT-directions, Xd 7^ 0, however, changes 9 and thereby 
decreases the interaction strength considerably. For a de- 
creasing transversal confinement we therefore expect the 
stable configuration to collapse by a displacement of the 
atoms in opposite A-directions in the polarized case. 

Tilted moments: The reasoning above is based on the 
assumption that the atoms are polarized, meaning that 
their dipole moments point in the direction of the electric 
field independent of the position of the atom. This is an 
oversimplification in case of Rydberg atoms in a strongly 
confining magnetic field configuration. In Section |II C| 
it has been shown that the application of a moderate 
electric field induces a dipole moment dp that is per- 
pendicular to the local direction of the magnetic field. 
In the Y-Z-plane dp points in the A-direction (we ne- 
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gleet here marginal dependences on Q). For finite X, by 
contrast, it has a finite Z-component, see App. |A 2| We 
thus concentrate once more on displacements in X and 
set Ys = Xd = in the following. The non-polarized 
case is illustrated in Fig. [4j 

In case of a displacement of the atoms in opposite di- 
rections, Xd > 0, the dipole moments of the two atoms 
include different angles with their connecting axis. They 
differ from the angles in the fully polarized case, 0p, by 
the additional tilt due to the local magnetic field direc- 
tion, ±A0. The interaction energy no longer depends 
on 1 - cos^iOp) but on 1 - cos(6'p + A6')cos(6'p - M) 
which is smaller than the former for Q <0p± Ad < tt/2. 
For electric dipole moments perpendicular to the local 
magnetic field axis, as considered here, the reduction of 
the dipole-dipole interaction energy for displacements in 
X]j is therefore smaller than in the fully polarized case 
discussed above. The curvature in that direction is thus 
expected to be still positive for shallower transversal con- 
finements. 

Symmetric displacements of the atoms in the same X 
direction, Xg ^ 0, do not change the dipole-dipole inter- 
action energy in the fully polarized case as stated above. 
For tilted dipole moments, however, the energy is reduced 
since the moments are no longer perpendicular to their 
connecting axis. This can be seen from the illustration in 
Fig. [4] The additional angle due to the orientation of the 
moments perpendicular to the local magnetic field axis 
(for Q = 0) explicitly reads 



G 

lS.0 ~ arctan ( —X 
B 



(41) 



The effects due to AO described above are therefore weak 
for typical parameter sets since for typical parameters the 
ratio G/B is small. A strong effect on the dipole-dipole 
interaction energy is expected for large ratios G/B. In 
this case, however, the transversal confinement is strong 
and the curvature in X-direction is positive on the Z-axis 
nonetheless. 

In order to verify the predictions above on how the sta- 
bility of the equilibrium configuration of the atoms gets 
lost, and the predictions for the configuration the atoms 
take when they step out of the Z-axis, we minimize the 
total two-atom potential for the coordinates Xs^ Xo, Ys 
and Yd for fixed symmetric displacements of the atoms 
in Z, Zd > and Zs — 0. As discussed above we expect 
the atoms to align on the Z-axis as long as the transver- 
sal confinement dominates the interaction. When the 
longitudinal confinement increases and/or the transver- 
sal confinement decreases, the atoms are expected to step 
out of the Z-axis in different X-directions, i.e., Xd ^ 0, 
Xs = Yd = Ys = 0, since the negative transversal gradi- 
ent at the equilibrium position due to the dipole-dipole 
interaction is largest in the Xu-direction. This change 
of the atoms' configuration is depicted in Fig. [s] where 
we show the computed Xd positions that yield minimal 
energy for fixed distances Zd as black dots into the two- 
dimensional section through the two-atom potential for 



different parameter sets (left plots). The series of plots 
in Fig. [5] shows the loss of the local minimum position for 
relaxing transversal confinement due to decreasing field 
gradients G (from top to bottom) . The vertical lines indi- 
cate the position of the equilibrium configuration on the 
Z-axis. The bar graphs on the right hand side show the 
minimal potential energy of the two atoms depending on 
the atomic distance in Z-direction, Zd (we still assume 
symmetric displacements in Z, i.e., Z5 = 0). As long as 
the equilibrium position on the Z-axis is a local potential 
minimum, the minimal energy concordantly exhibits an 
energetic barrier towards smaller distances of the atoms. 
The peak of this energetic barrier is located at one of the 
two saddle points of the potential that are located sym- 
metrically with respect to the longitudinal axis. When 
the local potential minimum is lost, these saddle points 
collapse into a single saddle point at the equilibrium po- 
sition on the Z-axis where the local minimum simultane- 
ously vanishes. This can be seen in the series of plots on 
the left hand side. The energy plots on the right hand 
side show that the energetic barrier also simultaneously 
vanishes. 

We find the minimal values of the potential energy at 
positions for which Xs is two orders of magnitude smaller 
than Xd, but nonzero. This is due to the permanent elec- 
tric dipole moments for vanishing electric field, which 
in turn are a signature of the finite size of Rydberg atoms. 
The origin of these moments is described in App. A 1 
and their properties in the X-Z-plane are depicted in the 
left subfigure of Fig. |4] The illustration shows that they 
point in the same direction as the electric dipole moments 
induced by the electric field {dp, middle plot) for positive 
displacements of the atoms in X. They are anti-parallel 
for negative X and they vanish on the Z-axis. Since d^ 
is significantly smaller than dp for all considered electric 
fields in this chapter, dj can be considered a correction 
to dp and their sum d is parallel to dp but larger or 
smaller in magnitude than dp for positive or negative X, 
respectively. This introduces an asymmetry in Xs and 
Xd into the dipole-dipole interaction energy and is hence 
responsible for the nonzero values of Xs for the position 
of minimal potential energy. This asymmetry is also re- 
sponsible for the negative values of Xd for the positions 
of minimal potential energy shown in Figs, [sj 

The stability of the equilibrium configuration is insen- 
sitive to changes in the electric field strength. The same 
is true for the transversal part of the center of mass wave 
function of the atoms. We also note that changing the 
principal quantum number n has no considerable effect 
on the stability of the equilibrium position as long as the 
requirements involving n can be met. The stability of 
the equilibrium configuration of the two atoms hence de- 
pends essentially on the magnetic field parameters B, G 
and Q but not significantly on F and n. 

In order to find an analytical stability condition involv- 
ing the magnetic field parameters we examine the curva- 
ture of the potential in Xpi-direction at the equilibrium 
position. For small ratios G/B (which is the case when 
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FIG. 5. (color online) Loss of the local minimum for decreasing magnetic field gradient G of the loffe-Pritchard trap, (a-c) 
Two-dimensional sections of the six-dimensional two-atom potential through the Xd-Zd plane, Xs ~ Ys = Yd = Zs = Q- 
The plot range of all three contour plots is ±150 kHz. (d-f) Minimal energy of the two-atom potential against Zd- Each point 
is computed minimizing the total potential for fixed Zd with the other center of mass coordinates as parameters. The Xd 
positions are shown as black dots in the sections on the left. Parameters: _B = lOG, Q = 6 x 1Q~^^B, F = 10"^" a.u., n = 30, 
G = 3Tm"^ [first row, (a,d)], G = 2Tm-^ [second row, (b,e)], G = l.STm"^ [third row, (c,f)]. 



the system is close to collapse), and assuming F ^ nG, 
we find the approximate expression for the curvature of 
the dipole-dipole interaction energy in X^j-dircction 



92 



dxl 



dd 



(42) 



which strongly depends on the distance of the atoms. 
Inserting the equilibrium distance ^_D,min from Eq. (31 1, 



where the dipole-dipole repulsion and the longitudinal 
confinement add to zero, we find 



dxl 



dd 



—6nQ 



(43) 



equilibrium 



For a stable configuration this anti-confinement has to 
be weaker than the transversal confinement due to the 
magnetic field at the equilibrium position, yielding the 
condition 



> UQ 



(44) 



As soon as the right-hand side of this inequality becomes 
as large as the left-hand side the two saddle points at 
the potential barriers join on the Z-axis and the local 
minimum is lost. Inserting the geometry parameter for 
the milhmeter trap in Ref. [3T], Q = 6 x 10" the 
stability condition takes the explicit form G[Tm"^] > 
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0. 17 S [Gauss]. 

The previously derived condition ( 40 1 is less restrictive 
than p4|) for all presented examples. This means that 



the collapse is interaction-induced and it is not due to 
the loss of magnetic confinement. 

For the existence of a stable configuration of the two 
atoms it does, in fact, not suffice to have a local mini- 
mum on the Z-axis, i.e., to meet the requirement (44 1. 



The potential around the local minimum must addition- 
ally be deep enough to accommodate at least the two- 
atom center of mass ground state. The height of the 
deciding potential barrier can be increased by tighten- 
ing the transversal confinement. For the exemplary pa- 



rameter set i? = 1 G, G = 1 Tm 



F = 10- 



a.u. 



Q = 6 X lO^^^B, the center of mass ground state en- 
ergy is Eq = 13.9 kHz, and the barrier height is 2140 kHz 
= 51.4 /zK. This temperature can typically be reached in 
a magneto-optical trap without further cooling. 



IV. EXCITATION SCHEMES 

The techniques that have been suggested and used to 
excite atoms into circular Rydbcrg states range from the 
microwave transfer method of Hulet and Kleppner [35], 
the crossed-fields method proposed by Delande and Gay 
|36j . and the RF field method proposed by Molander et 
al. [37j. The concern of this section, however, is not the 
excitation of single atoms into circular states but the ex- 
citation of two atoms in a loffe-Pritchard trap with an 
additional electric field directly into the equilibrium con- 
figuration that is stabilized by the dipolar interaction of 
the atoms. Solutions to this complication include (i) ex- 
ternally forcing the excitation to happen at the desired 
positions only, or (ii) changing the single-atom potential 
such that its minima coincide with the two-atom poten- 
tial minimum, or even (iii) adiabatically transferring the 
system to the desired state by varying the detuning of the 
laser during the excitation process. A possible implemen- 
tation of each of the options is described in the following. 
In principle all schemes are extendable to more than two 
atoms. 

The only coordinates that have to be externally im- 
posed are the Z-coordinates of the atoms since zero 
transversal displacement minimizes the energy for both 
the single-atom as well as the two-atom potential. The 
cleanest way to do this is to trap individual ground 
state atoms in two optical dipolc traps at the desired 
Z-positions. The trapping volume of such optical tweez- 
ers [38j can be made small enough (less than a /im in 
diameter) that only one atom can be captured in each 
trap [39]. These two atoms can then be excited using 
one of the methods described in Ref . [3S] . 

Another way of forcing the Rydberg atoms to be pro- 
duced at the desired Z-positions is to excite them from a 
cold ground state atomic cloud by two laser beams per- 
pendicular to the Z-axis that are focused next to each 
other to the desired equilibrium positions of the Ryd- 



berg atoms. This is possible if the equilibrium distance 
is considerably larger than the waist of the focused laser 
beams, which can be as small as one /im. Due to the 
strong Rydberg-Rydberg interaction, which yields an en- 
ergy shift within the excitation volume that is larger than 
the linewidth of the laser, only one atom can be excited 
within one of the laser beams. As the excitation can be 
located at any of the atoms in that region, however, the 
ensemble of atoms is excited collectively into a superpo- 
sition state called superatom |40J. 

The second solution involves the modification of the 
single Rydberg atom potential. This can be done by 
adding an extra wire on the X-axis to the Z-trap on a 
chip or, correspondingly, by adding an extra coil between 
the coils of a macroscopic loffc-Pritchard trap that are 
responsible for the loffe field. Both setups yield a dou- 
ble well potential with a variable barrier height and a 
variable distance of the potential minima. For vanishing 
electric field two Rydberg atoms can be excited indepen- 
dently from each other, one in the bottom of each well, by 
tuning the laser just under the energy of the minimum. 
In order to keep heating as low as possible, the magnetic 
barrier can now be substituted by the dipolar repulsion 
between the atoms by decreasing the current through the 
extra coil or wire and simultaneously increasing the elec- 
tric field strength. 

The circularization of the Rydberg atoms with a mod- 
ified adiabatic rapid passage method |41j . for example, 
can be completed within 5 ^s. The timescale of changing 
the magnetic field strongly depends on the configuration. 
If it is small enough, the described excitation scheme is 
scalable to produce more than two excitations, i.e., a Ry- 
dberg atom chain. This can be done by applying a mag- 
netic field gradient in Z-direction which tilts the trap 
and moves the stable Rydberg atom pair in Z-direction. 
The magnetic barrier can be ramped up again as to con- 
fine the pair in one of the wells. At the minimum of the 
other well an additional circular Rydberg atom can be 
produced. The two atoms in the first well mutually tune 
themselves out of resonance of the exciting laser due to 
their interaction. By ramping the magnetic barrier up 
and down, and exciting a Rydberg atom in the empty 
well every time as described, a stable chain of atoms can 
be produced in the trap. The procedure is restricted 
by the timescales of excitation, magnetic field switching 
and by the lifetime of the circular Rydberg atoms which 
is about 2 ms for n = 30 in the field-free case and scales 
with rt' . 

With both schemes mentioned above, the atoms are ex- 
cited into single-atom potential minima whose positions 
have to match the minimum of the two-atom potential 
that includes the interaction. Instead of artificially creat- 
ing single-atom potential minima outside the origin, one 
can adiabatically transfer a cloud of ground state atoms 
from the state with no excitation via the state with one 
excited atom at the origin to the stable equilibrium state 
for two atoms relying on the structuring effect of the 
dipole-dipole repulsion. The electric field must thus be 
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switched on at the beginning of the procedure. 

The idea is based on the dynamical crystalhzation ap- 
proach of Pohl et al. The starting point is a cold 
gas of ground state atoms that can be modeled as con- 
sisting of two-level systems. The two levels considered 
in Ref. |42j are the ground and the low angular momen- 
tum nSi/2 Rydberg state. Here, of course, we need the 
excited level to be the circular Rydberg state, which can 
be achieved by utilizing the rf-optical excitation tech- 
nique described by Cheng et al. [13] . Within the dynam- 
ical crystallization scheme, the coupling laser is detuned 
against the two-photon resonance. For large negative 
detunings the many-body ground state in the rotating 
frame of reference coincides with the initial state where 
all atoms are in the ground state. Increasing the detun- 
ings to positive values effectively lowers the energy levels 
of many-body states with one and two and more exci- 
tations. They cross at critical detunings Aj, A2, 
and states with 1, 2 and more Rydberg atoms are pop- 
ulated. The detuning is hence a control parameter that 
decreases the energy difference of adjacent number states 
|0), |1), |2),. . . with zero, one, two. . . excitations, respec- 
tively [44j . Since the laser couples the different number 
states, their energies undergo avoided crossings of sepa- 
rations 5\,. . . at the critical detunings. An adiabatic 
preparation of the states |1), |2),... is possible as long 
as the time in which the detuning of the laser changes is 
large compared to 1/(5°, l/5\,. . . . 

At the first crossing the initial state |0), with all atoms 
in the ground state, is directly coupled to the first excited 
state |1), with one Rydberg atom at the origin, which 
yields 5i ^ Vl. From |1) to |2), however, there is no 
direct laser coupling since the energetically lowest state 
with two Rydberg atoms, |2), is the stable equilibrium 
configuration described in Sec. |IIIB| with two Rydberg 
atoms symmetrically displaced from the origin. To go 
from |1) to 1 2), two off- resonant intermediate steps are 
required. First the central atom is de-excited and subse- 
quently the two Rydberg atoms at their equilibrium posi- 
tion are excited. A three-photon process is hence needed 
to come from |1) to |2). Assuming that |1) and |2) are 
resonant at time t, then the intermediate states are de- 
tuned by Aili. For the parameters used in Section IIIB 
(e.g., in Fig.jsl B = 30G, G = lOTm"!, Q = 6x10"^^ 
F = 2 X 10 a.u., n = 30) all detunings Afl^ are of the 
order of ~ 100 kHz. If the Rabi frequency is much smaller 
than the intermediate state detunings AQi, then the in- 
termediate states that couple |1) and |2) act as virtual 
levels for a resonant multi-photon transition. For larger 
Rabi frequencies, however, > Afi^, power broad- 

ening exceeds the intermediate state detunings and the 
states are coupled by consecutive one photon transitions. 



V. BRIEF SUMMARY 

In the present work we investigated the controlled trap- 
ping of two individual Rydberg atoms by means of a 



magnetic loffe-Pritchard trap that is superimposed by a 
constant electric field. The single-atom adiabatic poten- 
tials for such a field configuration have been derived and 
discussed. Including the interaction of the two Rydberg 
atoms, analytic expressions for the equilibrium positions 
of the two involved Rydberg atoms could be derived in 
the regime of a strong transversal confinement. As an in- 
teresting result, it turned out that the distance between 
the two atoms can be easily tuned without altering the 
involved trap frequencies by changing the applied elec- 
tric field. Loosening the restriction of a strong transver- 
sal confinement, on the other hand, leads to truly three- 
dimensional potential surfaces that, in principle, allow for 
the collapse of the system. The regime of stable trapping 
has been identified and the resulting adiabatic potentials 
were discussed. Possible routes to experimentally realize 
the proposed system have been outlined. 
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Appendix A: Computing electric dipole moment 
expectation values perturbatively 

In order to be able to base the perturbative treatment 
in the uppermost surface on the circular Rydberg state, 
the perturbation operator has to be rotated into the local 
direction of the magnetic field. 



Permanent electric dipole moments as finite size 
effect 



In the analytically diagonalizable case of high loffe field 
the electronic state corresponding to the uppermost elec- 
tronic adiabatic energy surface is the circular state with 
respect to the local field direction as quantization axis. 
For a finite gradient G the purely electronic finite-size 
term admixes states to the circular state that have 
opposite parity. The total wave function thus looses its 
definite parity and the matrix elements of the odd dipole 
operator er no longer vanish identically due to symmetry. 
This can result in a permanent electric dipole moment. 

Because the unperturbed state vector is analytically 
given in the rotated system, i.e., with the local mag- 
netic field axis as the quantization axis, the perturba- 
tion operator has to be rotated into this local frame, too. 
It is convenient to additionally replace the momentum 
operator in by angular momentum operators. This 
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can be done exploiting the energetic degeneracy of an n- 
manifold in a field- free environment. The commutator 
[xyz, Ha\ vanishes within an n-manifold, it is therefore 
{Lp\xyp z\'4>) — \{'p\{xLx — yLy)\ip), and the perturbation 
operator in the transformed frame becomes 

W = Un^U'' UX{xL^ - yLy + 3xS^ - SySy)U'f 

=: X c^j{R, B, G, Q) n{L, + 35^) , (Al) 

where A = G/S a.u. and TZa(} is the rotation associated 
with the transformation U , cf. Eq. (|9|. We write 7?."^ in 



Eq. ( Al I instead of TZap to recall the fact that the com- 
ponents of a vector operator transform in the rotation TZ 
like those of a vector in the rotation TZr"^ [45j. Both the 
coordinate vector r and the angular momentum opera- 
tors L and S are vector operators. 

We note that the perturbation operator, which is of 
purely electronic nature in the laboratory frame, now de- 
pends on the center of mass coordinates R through the 
coefficients Cij{R, B, G, Q) in the rotated frame. 

The first order energy correction to the uppermost cir- 
cular state IV'c) vanishes, 

Aei = (Vc|Ac,yr,(L, + 3Sj)\iPc) = . (A2) 

The first order correction in the wave function reads 

(1) _ {^p\Xcij ri{Lj + 3Sj)\(pi' 



\XipY 



(A3) 



Here we introduced the abbreviations Epq — E^ — E^ and 
Wpq — {ipp\Xcijri{Lj + 3Sj)\ipq), and we use the symbols 
\ip) = \n,l,m,ms) for the unperturbed hydrogenic elec- 
tronic states in energetic order, starting with the circular 
state constituting the uppermost surface. 
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2,n- 


2,-1/2) 



(A4) 

The states {|(/32), |<^3)}, the states {\lp4) , . . . ,\ipj)} and 
the states {jiySs), ■ • ■ , I'/'is)} are energetically degenerate 
in the limit B/G — > oo. The quantum numbers are given 
with respect to the local quantization axis which is the 
direction of the magnetic field. 



In order to co mpu te the matrix elements Wpi = Wp^^^ 
we rewrite the angular momentum 



defined in Eq. (A3 



operators with ladder operators, 

L.lVc) = ^(£++L_)|^,) = ^L^IV^c), 



(A5) 
(A6) 



where L^\l,m) — + 1) — m{m — 1) \l,m — 1), {fi — 
1 atomic unit). The only non- vanishing matrix elements 
in ( A3 1 are 



(^^|c„r,L,|Vc) = -\/2^^(^^|c„r,|^2) , (A7) 



(^^|c,^r,L,|Vc) = i-V2^^{'^lhyr,\ip2) , (A8) 
{^PslctzriL^lipcip) = (n- l)(</?3|ciz?'.i|V'c) (A9) 

due to the dipole selection rules Al = ±1, Am; = 0, ±1, 
and AttIs — 0. We proceed similarly with the spin oper- 
ators, 

S^lrus = ±i) = i |m, = T^) and (AlO) 
Sy\m, ^ = ±'- \ms ^ ■ (All) 

The only non-vanishing matrix elements involving the 
spin operators are 

{'PislcixriS^l'ipc} = ^{fislcixrilfj) , (A12) 

ifishyr^Sylipc) = ^{'Pi3\ciyr^\(p7) , (A13) 

ifslcizriSzlipc) = ^{fslczzrilipc) ■ (A14) 

Considering the following relations between the dipole 
matrix elements, 

{l\ m'\y\l, m) — ±j (/', m'|x|Z, m) for m' —m =p 1 
{I' , m' \y\l , m) — — {I' , m' \x\l , m) ior m' —m 
{l',m'\z\l,m) Sm,Tn' , (A15) 



we find the first order correction to the wave function, 
l^'Pi'^) = + icyz)^-^ ({n - ^)x3i + ^\/2n - 3z32^ 



1^) , 

Ei5 



1 



~l~ jy, {^xx ^yy ~^ ^''^xy / 2 



1^13; 

£^13, 5 



c„„ -I- 2iCxy)-^/2rl-^x<i2 
3 



"i~ \^xx Cyy 

+ 2icxy)-xi3j )-X . (A16) 



Here we introduced the notation Xij — {ipi\x\Lpj), yij — 
{LPi\y\ipj) and Zij = {(pi\z\ipj). The following explicit ex- 
pressions for the matrix elements can be deduced from 
the formulas for the radial and angular integrals involving 
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hydrogenic wave functions in 
3 



Z32 — 

Xl3 = 

X52 = 

2^13,7 = 



2 ' 
3 

3 

2^/2' 



3 /(n- l/2)(n-3/2) 



Vn- 2 



(A17) 
(A18) 
(A19) 

(A20) 



Note that the correction |A(/3^^'') to the circular wave func- 
tion \ipi) — lipc) has definite parity since {ifz), \ip5) and 
\1p13) have the same I quantum number. It is opposite 
to the parity of \tpc), however. The involved coefficients 
Cij{R, B,G,Q), defined in Eq. (Al), come from the in- 



verse rotation of r, L and S with TZap- Expressed via the 
rotation angles a and /3, the coefficients read explicitly 



= cos^ /? , 

= sin a sin /3 cos /3 , 
: — cos a sin /3 cos /? , 
: sin a sin /3 cos /3 , 
= (sin^ a sin^ /3 — cos^ a) , 
= — sin Of cos a(l -I- sin^ /3) , 
- — cos a sin /3 cos /3 , 
: — sin a cos a{l + sin^ /3) , 
:(cos^ a sin^ /3 - sin^ a) , (A21) 

where Cij = Cji. This is a general expression for the per- 
turbation operator (Al). The particular magnetic field 
configuration only enters via the explicit expressions for 
the angles a and f3. On the Z-axis {a = /3 = 0, U = 1) 



Cyx 
^xy 

Czy 

Cyz 
^zz 

Cji. 



all the coefficients but Cx 



1 and Cyy 



-1 vanish. The 



correction to the circular state in first order reduces to 

-|<^5) , 3 



|A^W)(0) 



>3 



51 



13,5 



(A22) 



for large n. The electric dipole moment expectation value 
therefore vanishes at the origin due to the dipole selection 
rules as will be described in the following. 

The electric dipole moment of the electronic state to 
second order in perturbation theory is found by comput- 
ing the expectation value of the dipole operator rotated 
into the local direction of the magnetic field, Wrll, for 
the perturbed state in the rotated frame, \ipi + A^?^^^), 

^{ipi + \^'~P\U^rU\ipi + X^'^P) 
^{ipi\U^rU\vi) + {X^^^\u^rU\X^^^^) 

+ 2 Re(((^i|[/V[/|A(^i^')) . (A23) 



The second line in Eq. (A23) vanishes due to the defi- 
nite parity of and \xlp^) and the matrix element 



in the third term simplifies due to the dipole selection 
rules (since Am = 2 for {ifi\UrU^\ipc,) and Atos = 1 for 
((pi|[/rC/t|(pi3)), 



-ia;i3 I 

V 



and it is 
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2^13 \ 




' Cx\ 


= 2 Re 










C-yz 
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/ 







where x :== 9n'^i2n^ - 3n - ^ ^•n? - lOn + 6 + 1) /(8Ai;) 
and A — G/Z. We can find the explicit spatial depen- 
dence of the electric dipole moment by expressing the 
angles a and /3 with the magnetic field components -Bi, 
to find the expression (13 1. For a loffe-Pritchard mag- 
netic field configuration, d-^ vanishes on the 2'-axis be- 
cause the magnetic field components and By are zero 
there. It is convenient for the symmetry analysis to also 
write down the explicit form of d~f for (5 = 0, 



d^(g = o) 




B^/G^) 
B^G^) 
(B/G) 



(A25) 



As can be deduced from Eq. (A24), the electric dipole 
moment expectation value is perpendicular to the local 
direction of the magnetic field. 



d^ B 



n 



^xz 
Ci, 



yz 

V 



B = Q 



(A26) 



2. Non-parallel moments in an electric field 

The Hamiltonian for the additional external electric 
field is Hp ^ qcf) = {xF^ +yFy + zFz) since F — — grad (p 
and q = —e {— —I in atomic units). According to the 
considerations in the preceding chapter the perturbation 
operator therefore reads 



Wf = U{r ■ F)U^ = UrU^ ■ F 



^F, ij Tj 



(A27) 



where the small parameter in the operator Wp is the 
modulus of the electric field, Xp = \F\. Consider- 
ing the Zeeman term dependence, fiB ^ \B\, the per- 
turbation parameter is the ratio of the field strengths, 
Xp = \F\/\B\. 

The first order energy correction due to Wp vanishes 
due to the dipole selection rules since Al — 0, 



Aj^efi^^) = (V'c|cf,.jF,r,#c) = 



(A28) 
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The first order correction to the circular state is 



{(Pp\cF,ijFirj\ip 



Ip 



Ei3 



fa) 



(A29) 



where we introduced Wp.pq = {^p\cF.ijFirj\ipq) and 
Epq = — E^. For an electric field in arbitrary di- 
rection the numerator in ( A29 ) reads 



Wf,31 = eX3i[Fx{cF^xx + iCF,xy) 

+ FyiCF,yy + F^ {cF,zx + iCF,zy)] (A30) 

and if we restrict our consideration to an electric field 
pointing along the X-axis we find 



IAfI^^^"^"^) = eXi3Fx{cF,xx +iCF.xy) 



1^3) 

Ei3 



(A31) 



where 

C-F,xx — cos /5 and CF.xy — sin o; sin ^. 

The second order energy correction due to the external 
electric field reads 



^2^,(2,F.) ^p^i^^^ ^ \Fv''i'''''\cF,xjr,Wi + 
=2Re [Fx{ipi\cF,xjrj\XFip^i''^"^) 

=2-=; Xi3{Cf,xx + '^F,xy) 



q f2 

„2,'„ iN/„„„«2 I „; 2, 

'4E 



13 



-n {n- l)(cos/3 + sin a sin ;9 ) 



(A32) 



For vanishing Q and with the approximate expression for 
the energetic separation between the coupling surfaces, 
Ei'i A£' « |B|/2, this reads 



A|e(^.-^) « ^F.V(n-l) |;+g^;^^,^, . (A33) 



The perturbative contribution to the uppermost surface 
due to an external electric field (F^,0,0) is thus positive 
and it is maximal on the z axis. 

Both the unperturbed wave function \ipi) as well as 

the perturbation |AfV?i ' ^ ) have definite parity. The 
electric dipole moment expectation value in the upper- 
most electronic energy surface is therefore, analogously 



to (A23I, 



dF = 2Re [{Lpi\U'^rU\XFipt' 



=7^;^(2x3l-^Re(c^^^ 
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I + sin a sin /3 
y sin /? 



9 Fx 



^13 



— 1) 
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sin a sin /? 
cos a ) 
y — sin a cos f3 ^ 

(A34) 



— BxBy 

—BxBz 

CF,y = sinasin/?, X13 = 



where we use cf,x — cos/3 

and the relation (A15|. In contrast to d^, dF 



does not vanish on the Z-axis but points in the direction 
of the electric field with dF^x = ^x\'3Fx/\B\. This is not 
true away from the Z-axis. The dipole moment does not 
point in the electric field direction there but stays rather 
perpendicular to the local direction of the magnetic field 
[alike d^, Eq. (A26)]. This can be deduced from the 



directional dependence in (A34) for arbitrary magnetic 
field configurations. 



dp B 



— BxBy 

\ —BxBz 



B = 



(A35) 



Addition of perturbatively calculated dipole 
moments 



For a non-zero external electric field the perturbation 
operator is Wtot = UHyW + UHfU'^ and the total 
first order correction to the wave function hence reads 
|A1) — |A^l)-y -I- |Af1)f- The calculation of the expec- 
tation value of the observable O in the perturbed state 
l^/jc + Al) yields accordingly 



(^e + AllOIV-c 

= 0.. 



Al) 
f Of 



2Re((A^l^|d|AFlF)) • (A36) 



The perturbations |A^1)^ and |Af1)f, Eqs. (AI6I and 
( |A29 1, involve the states It/^a), ly^s), |</5i3) which do not 
differ in their angular momentum quantum number /. 
They have the same definite parity since the parity of 
the spherical harmonics does not depend on the quantum 
number m. The mixed matrix element (A^l^|d|AFlF) 
therefore vanishes (the dipole operator is an odd opera- 
tor) and 



d= {Tpc + Al|r|-0c + Al) ^d^ + dF 



(A37) 



18 



In other words, within first order perturbation theory, and dp is equivalent to calculating the expectation 

adding the different dipole moment expectation values value for the combined perturbed state vector. 
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